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Surface Flows From Magnetograms 
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Abstract. Estimates of velocities from time series of photospheric and/or 
chromospheric vector magnetograms can be used to determine fluxes of mag- 
netic energy (the Poynting flux) and helicity across the magnetogram layer, and 
to provide time-dependent boundary conditions for data-driven simulations of 
the solar atmosphere above this layer. Velocity components perpendicular to 
the magnetic field are necessary both to compute these transport rates and to 
derive model boundary conditions. Here, we discuss some possible approaches 
to estimating perpendicular flows from magnetograms. Since Doppler shifts 
contain contributions from flows parallel to the magnetic field, perpendicular 
velocities are not generally recoverable from Doppler shifts alone. The induc- 
tion equation's vertical component relates evolution in B z to the perpendicular 
flow field, but has a finite null space, meaning some "null" flows, e.g., motions 
along contours of normal field, do not affect B z . Consequently, additional infor- 
mation is required to accurately specify the perpendicular flow field. Tracking 
methods, which analyze dtB z in a neighborhood, have a long heritage, but other 
approaches have recently been developed. In a recent paper, several such tech- 
niques were tested using synthetic magnetograms from MHD simulations. Here, 
we use the same test data to characterize: 1) the ability of the induction equa- 
tion's normal component, by itself, to estimate flows; and 2) a tracking method's 
ability to recover flow components that are perpendicular to B and parallel to 
contours of B z . This work has been supported by NASA Heliophysics Theory 
grant NNG05G144G. 



1. Why study surface flows from magnetograms? 

The large length scales and relatively high conductivity of the plasma in the 
solar corona imply that, to a good approximation, magnetic flux is frozen to 
the plasma there. Consequently, the coronal magnetic field is "line-tied" to the 
plasma in lower atmospheric layers where hydrodynamic forces can be stronger 
than Lorentz forces — the photosphere and lower chromosphere — and coronal 
evolution is strongly coupled to evolution in these layers. Accordingly, observa- 
tions of magnetic field evolution below the Sun's corona — typically, sequences 
of photospheric or chromospheric magnetograms — provide crucial tools under- 
stand coronal evolution. 

Usually, vector magnetograms are more useful than line-of-sight (LOS) mag- 
netograms for studying coronal evolution, because information derived from LOS 
measurements alone will not, in general, be consistent with the actual magnetic 
field, which has field components both parallel and transverse to the LOS. Al- 
though time series of vector magnetograms have historically been rare, SOLIS 
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(jHennev et al. 20021 ). the Solar Optical Telescope (SOT, Tarbell 2006) on Hin- 
ode, an d the Solar Dynamics Observatory's Helioseismic and Magnetic Imager 
(HMI) (jScherrer and The HMI TEAM 20051 ). should dramatically improve pho- 
tospheric vector magnetogram spatial and temporal coverage in the near future. 

Several tec hniques have been developed to derive flows from time series of 
magnetograms (Chae 2001 ; Kusano et al. 20021 ; IWelsch et al. 20041 ; lLongcope 20041 ; 
Georgoulis and LaBonte 20061 : ISchuck 20061 1. Estimated flows at the base of the 



corona can be used t o derive the fluxes of magnetic helicity, energy, and free en- 
ergy into the corona, ( Demoulin and Berger 20031 ; iPariat et al. 20051 : 1 Welsch 2 006). 
Further, flow estimates can be used to provide time-dependent boundary condi- 
tions for data-driven simulations of coronal magnetic field evolution. 

This paper briefly reviews progress on estimating surface flows from mag- 
netogram sequences, and demonstrates some aspects of the problem with test 
data. 



2. Progress in estimating surface flows from magnetograms 

How have velocities been derived fro m vector magnetograms? Ch ae (2001) ap- 
plied local correlation tracking (LCT) ([November and Simon 19 88) to LOS pho- 
tospheric magnetograms to determine the proper motions of magnetic features 
on the magnetogram surface, and assumed the inferred flows tJ LCT ) were es- 
timates of the horizontal plasma velocities v^. (Here the h subscript denotes 
a vector's components tangential to the magnetogram surface. We avoid the 
subscript t, which has been used to refer to vector components both transverse 
to the LOS and tangential to the solar surface.) Assuming that the observed 
evolution of the photospheric magnetic field is governed by flows according to 
the ideal induction equation, 

<9B 

— = -c(VxE) = Vx(vxB) , (1) 

Kusano et al. (2002) proposed using the component of equation (pQ) normal to 
the magnetogram, 

8B 

= z • V x (v x B) = -V • (y h B z - v z B h ) , (2) 

to derive three-component velocity fields — v x ,v y ,v z . Here, we have defined 
the magnetogram surface to be the horizontal plane containing x and y, with 
a vertical normal z. The three components of v cannot completely determined 
from d2|) alone, so more data or assumptions are required to close the system for 
v. 

What about the horizontal components of equation ([I])? As Kusano et 
al. (2002) noted, only the normal component of equation ([1]) is completely spec- 
ified by vector magnetic field measurements from a single atmospheric layer; 
the other components of equation (pQ) contain vertical derivatives of horizontal 
magnetic field components, and therefore require measurements of the vector 
magnetic field at a different hei ght in the atmosphere ( e.g., the chromosph ere) 
which are only rarely available ( Leka and Metcalf 20031 ; [Metcalf et al. 20051 ) . 
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Kusano et al. (2002) assumed LCT velocities, u, to be equivalent to to 
close the system. Demoulin and Berger (2003) argued that tracked motions of 
magnetic flux on the solar photosphere, u, result from the combined effects of 
horizontal plasma velocities transporting vertical magnetic fields and vertical 
plasma velocities transporting horizontal magnetic fields, via 

uB z = w h B z - v z B h . (3) 

The distinction between apparent motions of flux and plasma velocities led 
Welsch (2006) to term u the "flux transport velocity." Demoulin and Berger 
(2003) suggested that LCT could be used to infer u, but not directly. We 
note that, in addition to ideal flu x transport, d iffusive effects can also lead to ap- 
parent flux transport velocities ( Welsch 20 06). although we ignore these effects 



in the present work. 

Since the seminal work of Kusano et al. (2002), still more techniques have 
been developed that determine velocities from vector magnetograms. Welsch 
et al. (2004) used equation (J3j) to combine Fourier LCT (FLCT) results with 
equation ([2]) to determine a photospheric flow field in a method they termed 
inductive local correlation tracking, or ILCT. Longcope (2004) developed the 
minimum energy fit (MEF), which finds the photospheric velocity field that is 
strictly consistent with equation ([2]) and that minimizes a penalty function, e.g., 
the integrated square of the three-component photospheric velocity. Georgoulis 
and LaBonte (2006) extended the minimum structure method of Georgoulis, 
LaBonte, and Metcalf (2004) to the problem of velocity determination, in a 
method they termed minimum-structure reconstruction (MSR). Schuck (2005) 
showed that, formally, LCT is not consistent with the induction equation's nor- 
mal component, which can be expressed as a continuity equation; instead, LCT 
is consistent with the advection equation. Building upon the "differential LCT" 
(DLCT) method developed by Lucas and Kanade (1981), Schuck (2006) de- 
veloped the differential affine velocity estimator (DAVE), which employs least- 
squares fitting to solve the continuity-form of equation ([2]) for u and its spatial 
derivatives. All of the methods listed above can be applied to chromospheric as 
well as photospheric magnetograms. 

Flow estimatio n techniques that explicitly em ploy equation ([2]) have been 
termed "inductive" ( Georgoulis and LaBonte 20061 ); examples include IM, ILCT, 



MEF, and MSR. Tracking methods, e.g., FLCT, DLCT, and DAVE, use the evo- 
lution of structure in B z to qu antify proper motions, and may also be referred to 
as "optical flow" techniques ( Schuck 20061 ). Tracking techniques usually apply 



a windowing function, centered on each pixel tracked, to derive the optical flow 
in the neighborhood of that pixel. Since this window imposes a scale length, 
tracking methods have bee n criticized for their selective s ensitivity to flows on 
the imposed window scale ( Georgoulis and LaBonte 20061 ). The tracking tech- 



niques mentioned thus far might be termed "Eulerian methods," as they estimate 
velocities over a set of pixels from an image pair. In contrast, "feature track- 
ing" techniques ( DeForest et al. 20 07). in which discrete "features" are identified 



from structure in B z maps and followed in time, might be termed "Lagrangian" 
techniques. We will not discuss feature tracking techniques further in this work. 

Generally, velocity estimation techniques are susceptible to errors in mag- 
netograms, because changes in the inferred magnetic field, dtB z , are assumed to 
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arise from flows via equation ([5]). Simply put, when noise or systematic errors 
introduce spurious fluctuations 5B Z in B z , spurious velocities are derived. In 
particular, in regions where \B Z \ is small, 5B Z due to noise can be on the order 
of changes AB Z due to actual field evolution. For this reason, tracking methods 
are expected to estimate flows more accurately in regions where \B Z \ is large 
(barring magnetograph saturation effects) than where \B Z \ is small. 

The recent proliferation of velocity estimation methods led Welsch et al. (2007) 
to test seven routines' ability to reproduce known flows. They extracted several 
pairs of "synthetic magnetograms" from ANMHD simulations of a rising flux 
rope in the upper solar convection zone, and used several methods methods - 
two LCT codes, DAVE, IM, ILCT, MEF, and MSR — to estimate flows respon- 
sible for magnetic evolution between each pair. They then compared properties 
of the estimated flows with those of the "true" flows from the MHD code. (These 
flow reconstructions were not blind: those deriving flow estimates had access to 
the true flow.) They found that MEF, DAVE, FLCT, IM, and ILCT performed 
similarly by many measures, but that MEF estimated the fluxes of magnetic 
energy and helicity quite well. The other methods tested estimated the fluxes 
of magnetic energy and helicity through the magnetogram layer poorly. 

While the MHD simulations of field evolution in the solar interior used 
by Welsch et al. (2007) did not accurately model photospheric field evolution, 
they did provide a valuable tool for testing velocity estimates. In the near 
future, blind tests, with simulated data that more accurately model photospheric 
field evolution, are planned, as are tests of the sensitivity of flow estimation 
methods to noise. (Application of these methods to synthetic chromospheric 
vector magnetograms will also be investigated.) 



3. Theoretical considerations 



As noted in the discussion of equation ([2]), dtB z alone cannot fully specify the 
three components of v. We note that both e yolution in B (equation [J) and the 



ideal, upward fluxes of magnetic energy S z ( Demoulin and Berger 20031 ) . 

S z = ^ J dA'z • [(v x B) x B] = ^- j dA' (uB z ) ■ B h , (4) 

and magnetic helicity (Demoulin and Berger 2003; see also Pariat et al. 2005) 

dH A /dt = 2 J dA' z • [A P x (v x B)] = -2 J dA' (A P ■ u)B z (5) 

are all linear in (v x B). (In equation [5], = [d x , d y , 0] T , z • x Ap = B z 
and Vfe ■ Ap = z ■ Ap = 0.) This dependence means that flows along the 
magnetic field, vn, do not lead to any evolution in B, and do not transport 
either magnetic energy or magnetic helicity across the magnetogram surface. 
Hence, the components of v perpendicular to B, 

v ± • B = , (6) 

are of practical interest for most applications of velocity estimation. This equa- 
tion, combined with equation ([2]), now provides two equations for the three 
unknown components of v. 
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In this section, we briefly consider possible remedies to the remaining un- 
derdetermination, as well as its implications. 



3.1. Can Doppler data determine surface magnetic flows directly? 

Doppler measurements of the velocity of the magnetized plasma, best deter- 
mined from Stokes V, Q, and/or U profiles (not J; see Chae et al. [2004]), can- 
not fully determine whether the inferred plasma flows lie either parallel to or 
perpendicular to the magnetic field if -Blos 0- This degeneracy is illustrated 
schematically in Figure [TJ 
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Figure 1. When -Blos 7^ 0, a flow vj_ perpendicular to B can produce the 
same LOS velocity as a flow vii parallel to B. This means LOS flows inferred 
from Doppler shifts are insensitive to the relative orientation of v and B. Flows 
parallel to B do not lead to evolution in B, and do not transport magnetic 
energy or helicity. The evolution of B z at the photosphere, however, can be 
used to estimate vj_ . The presence of a field component perpendicular to the 
page (neglected in this figure) would not affect this analysis. 



Doppler shifts do unambiguously determine the perpendicular flow when 
Blos = 0. Near disk center, this is approximately s atisfied along po larity 
i nversion li nes, where B z changes sign, as noted by ( Chae et al. 20 04) and 
( Lites 20051 ). Doppler observations of have shown that some polarity inversion 
lines can also be velocity inversion lines (Deng et al. 2006), which can be inter- 
preted as the manifestation of siphon flows moving along field lines that arch 
over the polarity inversion line. 

Doppler measurements can, however, be combined with an estimate of the 
perpendicular velocity vj_ from one of the techniques described above to recover 
the flow parallel to the magnetic field: subtracting the projection of vx onto the 
LOS from wlos gives the LOS component of v\\, which can be divided by the 
cosine of the angle between B and the LOS to give vn . Georgoulis and LaBonte 
(2006) have employed this approach in an observational study of active region 
flows, and Ravindra and Longcope (2007) have investigated the use of Doppler 
shifts in a theoretical study using the data from Welsch et al. (2007). 

In the following discussion, any reference to v should be assumed to refer 

to Vj_. 
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3.2. Inductive flows 

One can employ a Helmholtz decomposition of equation ([3]), 

uB z = \ h B z - v z B h = -VhX - V ft x , (7) 

to express uB z in terms of "inductive" and "electrostatic" potentials ( Longcope 20041 ) 
X and ip, respectively. (We note that, assuming the ideal Ohm's Law, these po- 
tentials can be expressed as sources of a horizontal electric field, 

B h = c-^z x uB z ) = c-\V h x X z) - c- l V h 4> .) (8) 

Inserting equation ([7]) into equation ([2]) yields a Poisson equation for Xi 

d t B z = V 2 hX , (9) 

meaning the evolution of B z between a pair of sequential magnetograms, AB Z /At, 
specifies x- This does not, however, constrain if) in any way. 

Equation ([U]) can be solved to determine the purely inductive flux transport 
rate, i.e., one that assumes ip = 0, 

UI B Z = -V hX , (10) 



where the subscript / denotes that the flux transport velocity is inductive. In 
principle, equation (fTO]) can be substituted into equation ([3]) and combined with 
equation ([U]) to determine to determine v. In practice, though, solving this 
system of equations for v requires division by B z , which weights regions of weak 
\B Z \ more strongly than regions of large \B Z \, and, as discussed in fjLj flows are 
expected to be poorly estimated in such regions. For some applications, however, 
an estimate of uB z is sufficient; for instance, equations © and © imply that 
uB z can be used to estimate the fluxes of magnetic energy and helicity into the 
corona without determining v. 

The ANMHD data used by Welsch et al. (2007) can be used to demonstrate 
the ability to estimate flow properties from the induction equation alone. In 
Figure [21 we show a snapshot of the magnetic (grayscale background for B z 
and white vectors for B/J and velocity (black/white contours for v z , and black 
arrows arrows for v/J fields in a cross section of the simulation domain near the 
top of the convection zone. We note that the strong fields present (~ 8000 G) 
are appropriate to the high-/? solar interior, but are substantially stronger than 
observed photospheric field st rengths. These dat a are described in much greater 
detail in Welsch et al. 2007 dWelsch et al. 20071 ). 

Using the AB Z between two synthetic magnetograms 125 seconds before 
and after the magnetogram shown in Figure [21 we solved equation ([9]) with a 
Fourier technique to compute uiB z . Figure [3] shows a scatter plot comparing 
ANMHD's flux transport rate, U iB z , with ujB z , along with the linear corre- 
lation coefficients for the x and y components. These correlation coefficients 
compare quite favorably with results of the tests by Welsch et al. (2007). We 
also computed a least-absolute-difference linear fit of |u/i?^| as a function of 
IU-RjI, and found a slope of 0.4, when 1.0 would be ideal; evidently, inductive 
flux displacements underestimate ANMHD's true flux displacements. 
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B & V from ANMHD Simulations 
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Figure 2. A snapshot of the magnetic and velocity fields in a cross section 
of an ANMHD simulation of a flux rope rising through the convection zone, 
extracted from near the top of the domain. The grayscale background shows 
B z , white vectors show B; l7 black (white) contours show smoothed downflows 
(upflows) of v z , and black arrows show v/j. Only components of v perpendic- 
ular to B are shown. 

In addition, we calculated the estimated fluxes of magnetic energy and 
magnetic helicity, separately trying the methods of Demoulin and Berger (2003) 
and Pariat et al. (2005) for the helicity flux. We found that ujB z determined 
by our Fourier method recovered 30% of the energy flux, but only 10% of the 
helicity flux. If, however, the estimated helicity flux was summed only over 
strong-field pixels — in which \B Z \ exceeded 5% of max(|i? 2 |) — the Demoulin 
and Berger (2003) approach recovered 20% of the helicity flux from all pixels, 
while the Pariat et al. (2005) approach recovered 15% of the helicity flux from 
all pixels. While these results compare favorably with most of the velocity 
estimation methods tested by Welsch et al. (2007) , they reveal that inductively 
determined flux transport rates lack essential information about the fluxes of 
energy and helicity — presumably related to tp being unconstrained by equation 
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Flux Transport Veloci;y frorr AB Z /At Alone 




Figure 3. A scatter plot comparing ANMHD's flux transport rate, UjB z , 
with that derived from equation ([9]), u/i? z , along with the linear correlation 
coefficients for the x and y components. The solid line is not a fit; it is the 
desired slope, and is shown to emphasize departures from that slope. 



The inductive potential x used in these test was determined with a Fourier 
technique, which assumes periodicity, and solves for \ on all pixels, including 
those with weak vertical field. We ran the same tests with an alternative poten- 
tial function, which satisfies homogeneous boundary conditions over a masked 
region that includes only pixels with |I? Z | > 170 G. Derived by B. Ravindra and 
D.W. Longcope, x' was used for the inductive component of MEF flows in the 
tests conducted by Welsch et al. (2007). We found that, compared to flux trans- 
port rates from x> flux transport rates derived from x' were more poorly corre- 
lated with ANMHD's flux transport rates, underestimated |U-B 2 | more severely, 
and recovered less of the vertical fluxes of magnetic energy and helicity. We note 
that MEF overestimated |U1? 2 | in the tests performed by Welsch et al. (2007) . 
This implies that the superior performance of the MEF approach at estimating 
fluxes of magnetic energy and helicity in Welsch et al. (2007) derives from MEF's 
specification of ijj. 
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3.3. Sensitivity of flow estimation to VhB z 

As Demoulin and Berger (2003) observed, the underdetermination of v by equa- 
tion ([2D implies that there exists a class of flows that cause no evolution in B z , 
but that can inject large amounts magnetic energy and helicity. Flows that do 
not alter B z satisfy 

Vh,oBz ~ v z>0 B h = u B z = -V /t x ipz . (11) 

Because these flows lie in the null space of equation ([2]), we call them "null 
flows," and denote them with a zero subscript. While flows along contours of 
B z — "contour flows" — are a well known subset of null flows, modelers have 
employed other null flows; see, e.g., ( Lynch et al. 20051 ). 



Taking the curl of equation (|7|) yields 

V h x uB z = • (12) 

A flux transport velocity, u^ EST \ estimated from any optical flow technique 
(e.g., DAVE or FLCT) can be u sed with (fT2l) to find ip, and equation © to find 
X- This is the essence of ILCT ( Welsch et al. 2004] ); but it goes a step further, 



and uses equations ([7]) and §5§ to determine v. Since, however, optical flow 
techniques ultimately depend on AB Z /At, and flows associated with tp do not 
alter B z , this procedure is expected to be insensitive to such flows. 

We can use FLCT flows to illustrate the difference in sensitivity of one 
velocity estimation technique to flows along contours of B z versus flows along 
gradients of B z . In Figure El we plot selected contours of ANMHD's average 
B z over a grayscale image of the change AB Z over the same 250 second time 
interval used to estimate ujB z in ^3.2.I We also plot ANMHD's instantaneous flux 
displacement, XJB Z , decomposed into components along contours (gradients) of 
of B z with white (black) vectors. Flow estimation techniques that depend upon 
dtB z are expected to be insensitive to flows that transport flux along contours 
of B z (white vectors). 

Using u from FLCT in equation (fT2j) to find ip, and equation © to com- 
bine x from equation (J9j) with tp., we attempted to reconstruct ANMHD's flux 
displacement, TJB Z , from Figured! In Figure we show scatter plots of flux 
displacements from ANMHD and our estimates of these (denoted TJB Z and uB z , 
respectively) along VhB z (horizontal dashes) and contours of B z (vertical line 
segments). Flows along V^-Bz are recovered more accurately than flows along 
contours of B z , as quantified by the linear correlation coefficients shown on the 
plot. The fact that contour flows are even partially recovered presumably re- 
sults from variations in the spatial structure of B z in the neighborhood of each 
tracked pixel. 



4. Discussion 

We have briefly reviewed central concepts regarding surface flow estimation from 
magnetograms. In addition, we have presented the results of simple tests which 
demonstrate that the induction equation's normal component, equation ([2]), can 
be used to quantitatively estimate flux transport rates. Fluxes of magnetic en- 
ergy and helicity derived from these estimates are, however, likely to possess 
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Figure 4. The change AB Z in the normal magnetic field is shown in 
grayscale, along with selected contours of B z (white for B z > 0, black for 
B z < 0). ANMHD's instantaneous flux displacement, XJB Z , decomposed into 
components along contours (white vectors) and gradients (black vectors) of 
B z , is overplotted. Flow estimation techniques that depend upon dtB z are 
insensitive to flows that transport flux along contours of B z (white vectors). 



significant systematic errors. We have also demonstrated that flow estimation 
techniques that depend upon evolution in B z alone are insensitive to flux trans- 
port along contours of B z , compared to flux transport along VhB z . Equations 
([3]) and ©, which depend on the dot products of uB z with B/j and with Ap 
respectively, combined with the sensitivity of uB z to VhB z , suggest that if 
and/or Ap lie primarily along Vh,B z , then the fluxes of magnetic energy and/or 
helicity can, in principle, be recovered accurately from AB Z /At. If, in contrast, 
B^ and/or Ap lie primarily along contours of B z , then the fluxes of magnetic 
energy and/or helicity probably cannot be recovered accurately from evolution 

in B z alone. 

We note that, as discussed in (jWelsch et al. 20071 ). the ANMHD data used 
in the tests presented here differ from actual magnetograms in significant ways, 
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Estimated Flux Transpor: on Contours &c Gradients 




ANMHD UB Z on Contours (I) & Gradients (-) [G cm s~ 1 ] 

Figure 5. Scatter plots of ANMHD's XJB Z and estimated uB z are shown. 
The vectors have been decomposed into components along gradients in B z 
(plotted with -)and contours of B z (plotted with |). Flux transport along 
gradients of B z is recovered more accurately than transport along contours 
of B z , as quantified by the linear correlation coefficients shown for each com- 
ponent. The solid line is not a fit; it is the desired slope, and is shown to 
emphasize departures from that slope. 

so the properties of flows estimated from actual magnetograms will probably 
differ substantially from the properties of flows estimated from ANMHD data. 

Acknowledgments. Data from MEF velocity estimates were courtesy of B. 
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